Stability of spontaneous quantum Hall state in the Triangular Kondo-lattice model 
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We study the behavior of the quarter-filled Kondo lattice model on a triangular lattice by combining a zero- 
temperature variational approach and finite-temperature Monte-Carlo simulations. For intermediate coupling 
between itinerant electrons and classical moments Sj , we find a thermodynamic phase transition into an exotic 
spin ordering with uniform scalar spin chirality and (Sj) = 0. The state exhibits spontaneous quantum Hall 
effect. We also study how its properties are affected by application of an external magnetic field. 

PACS numbers: 71.10.Fd, 71.27.+a, 73.43.-f 



The very broad spectrum of physical phases and responses 
of strongly correlated materials is rooted in the simple fact that 
electrons carry both charge and spin. The interplay between 
these elementary degrees of freedom gives rise to unconven- 
tional forms of superconductivity, multiferroic behavior, giant 
magnetoresistance, and heavy fermion physics. It also entails 
the possibility of strong magneto-electric effects. In insula- 
tors, the electric polarization can be induced or modified by 
applying a magnetic field, or by coupling to certain magnetic 
orderings fl]]. In metallic systems, the magneto-electric ef- 
fects are manifested in the conductivity tensor. The giant 
and colossal J3l magneto-resistance effects are examples in 
which the diagonal part of the conductivity tensor is dramat- 
ically modified. In addition, magnetic ordering can lead to 
changes in the off-diagonal components of the conductivity 
tensor. For instance, the "anomalous Hall effect" observed in 
the presence of ferromagnetic order, even in the absence of 
externally applied magnetic field, H = 0, has been a subject 
of active research over the past several decades 

Recently, a residual Hall effect was observed in the metal- 
lic pyrochlore system Pr2ir2C>7, in zero magnetic field and in 
the absence of uniform magnetization |8y. This observation 
strongly suggests that the effect is caused by a magnetic struc- 
ture with non-zero average scalar spin chirality, which for a 



single triangular plaquette is (x 



ijk ) = (Si ■ S, S; i|. 

Scalar spin chirality breaks time reversal and parity symme- 
tries and can be stabilized even in the absence of usual mag- 
netic ordering: (Si) = [10]. The symmetry properties of 
this order parameter lead to unusual magneto-electric effects 
both in metals |@, U - 15 ] and in insulators fTrl . The origin 
of the magneto-electric coupling lies in the Berry phase Il7ll 
that electrons accumulate as they traverse closed paths in the 
real space. The Berry phase is equal to half of the solid angle 
subtended by the electron spin as it moves around the path; 
e.g., in the continuum limit of smooth magnetic textures, it 
is proportional to the integral of the scalar spin chirality (or 
Berry curvature) over the area enclosed by the loop. For a 
given spin species (locally parallel or antiparallel to the mag- 
netic texture), the effect of the Berry phase is equivalent to 
an orbital magnetic field, and therefore can lead to finite Hall 
conductivity even in absence of an external magnetic field. 

Under the special circumstances of commensuration be- 
tween the strength of the effective magnetic field (i.e., Berry 



curvature of magnetic texture) and the itinerant electron den- 
sity in quasi-2D systems, there is an exciting possibility of ob- 
taining spontaneous quantum Hall effect (SQHE), i. e., quan- 
tized Hall response in the absence of magnetic field or uniform 
spin polarization. In a recent work 111 lh . we have shown that 
such insulating chiral phase is indeed stabilized at zero tem- 
perature, T = 0, in the weak-coupling regime of a triangular 
Kondo lattice model (KLM) for a 3/4 filled conduction band 
(p = 3/4). This result was based on the perfect nesting prop- 
erties of the non-interacting Fermi surface for p = 3/4, which 
leads to a "3Q ordering" of classical local moments at T = 0, 
with large uniform scalar spin chirality. 

In the present work we explore the stability of chiral order 
in a 2D Kondo lattice model (KLM) with respect to thermal 
fluctuations and magnetic field. It is well known that a con- 
tinuous symmetry cannot be broken at finite temperature in 
2D as long as the interactions are of short range 11811 . How- 
ever, a non-coplanar magnetic ordering breaks a discrete Z2 
symmetry, in addition to the continuous symmetry, that cor- 
responds to the two disconnected 5*0(3) sectors of 0(3] ro- 
tations needed to parametrize the order parameter Il9l I2O1I . 
The Ising component of this order parameter (chirality) can 
survive at finite temperature leading to a chiral spin-liquid: 
(Xijk) 7^ and (Sj) = 0. The interplay between the contin- 
uous and the Ising degrees of freedom in such chiral liquid is 
non-trivial and can lead to deviations from the Ising universal- 
ity class and a reduction of the transition temperature l2lLl22ll . 

To study this problem we combined a zero-temperature 
variational approach with finite-temperature Monte Carlo 
(MC) simulations of the KLM with classical local moments. 
We demonstrate that the 3Q phase with uniform chirality 
(UCP) of Ref. UM is stable not only at p — 3/4, but also 
at p = 1/4 in the intermediate coupling regime 1 < J/t < 4 
(t is the hopping amplitude for the conduction electrons and 
J is the exchange coupling to the localized moments). This 
result agrees with the small-unit cell variational calculations 
at T — presented in Ref. 112311 . Our MC results pro- 
vide evidence for a first order thermodynamic phase transition 
into the UCR The obtained ordering temperature at J = 2t, 
T c ~ 0.03t, is relatively small compared to J and t, likely 
because of the strong frustration and the suppression of the 
chiral order by the continuous fluctuations of magnetization. 
Nevertheless, the value of T c may be high enough for observ- 



ing SQHE near room temperature in transition metal oxides 
(assuming that t is of order leV 12411 ). Finally, we consider 
the effect of a magnetic field acting on the local moments both 
on the stability of the T = UCP phase and the spontaneous 
quantum Hall effect, as well as on the the value of T c . 

We consider the Kondo-lattice Hamiltonian on a triangular 
lattice with periodic boundary conditions, 



H 



= -*£ 



C la C 3<7 



where c^ a (cj a ) is the creation (annihilation) operator of an 
electron with spin a on site j, Sj is a classical Heisenberg 
spin with |Sj| = 1, cr^ = (af v , a^, a*) is a vector of Pauli 
matrices, and indicates that I and j are nearest-neighbor 
sites. Since the sign of J is irrelevant for classical moments, 
Sj, we will assume J > for concreteness. 

The state of interest is characterized by the local scalar spin 
chirality Xijk = &i ' Sj x Sk, and its Fourier transform: 
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FIG. 1: (color online) The 4-sublattice spin structures, (a) Triangu- 
lar lattice, with a — d indicating each sublattice (b) High-symmetry 
four-sublattice spin structures. The nomenclature follows the one for 
magnetic ordering on a tetrahedron with vertices (abed). 



The index a denotes each triangular plaquette, and N = L 2 
is the total number of lattice sites. The global order parameter 
for the UCP is (xo)- Local spin correlations are described by 
the spin structure factor 



^ k ) = ^E( s '- s ^ elk ' (r 

3,1 



(2) 



which is useful for characterization of the T = spin order- 
ing. We note here that for the perfectly ordered UCP, same 
as the "all-out" phase in Fig. Q] simple arithmetic shows that 
Xo ^ 0.59 and 5(0) = 0, while x q = and 5(0) = N for 
the fully polarized ferromagnetic state. 

The MC simulation samples the space of all possible classi- 
cal spin configurations. For each configuration {S}, the elec- 
tron eigenstates are found by exact diagonalization of H. To 
study the thermodynamic properties at a fixed filling factor p, 
we perform a Legendre transformation of the free energy in 
the grand canonical ensemble. The resulting free energy is 

F(N,/3) ee -~lnS, E = J2W({S}), 
1 {S} 



W({S}) 



(3) 



The chemical potential M{s} is adjusted for a given configura- 
tion {S} such that the total number of electrons, N e = J2\ fx, 
is the same for every spin configuration. Here, e\ is A- 

th eigenvalue of H({S}), fx = (e p { tx ~^ s i} + 1) , and 
P = 1/T. In this work we compute the energy density, 
e = N^ 1 dp(fiF), and the specific heat, c = dxe. Our nu- 
merical results are based on the single spin flip MC dynamics, 
with the flip decided by applying the Metropolis algorithm. 



We perform 10000 MC sweeps for each processor and esti- 
mate the statistical errors using 8 mean values typically. 

First we consider the T = case by means of a variational 
calculation in which we minimize the total energy over sev- 
eral highly symmetric four-sublattice spin structures shown in 
Fig. 03 12611 . The "all-out" and ferromagnetic configurations 
are the only two states that are stabilized as a function of J/ 1 
for p = 1/4. Figure |2{ a) includes a comparison between the 
corresponding energy densities e({S}) = L~ 2 J2\ e \f\ f° r 
p = 1/4. The solid lines are the results for L = 512. The 
all-out structure has the lowest energy for < J/t < 4.9, 
while the ferromagnetic configuration becomes the minimum 
energy state for J/t > 4.9. A similar comparison for L = 8 
(dashed lines) indicates that important size effects appear in 
the weak coupling regime, J/t < 1, for small values of N 
(see the inset of Fig. |3a).) The energy of the ferromagnetic 
structure is lower than the energy of the all-out structure for 
L = 8 and small values of J/t < 1, although this is not 
true in the thermodynamic limit l27ll . The all-out structure, 
which is the lowest energy variational state for J/t < 4.9, has 
a net uniform scalar spin chirality. Therefore, it has to pro- 
duce a spontaneous Hall effect. By explicit calculation of the 
Hall conductivity from the electron band structure, we find the 
quantized value a xy = ±e 2 /h for 0.7 < J/t < 4.9. 

The MC calculations at low temperatures allow us to test 
the results of the variational approach, not being limited to the 
2x2 magnetic unit cell size. Figure|2|b) shows our MC results 
for the J/t dependence of Xo an d S(0) at a very low temper- 
ature T/t = 10 -4 and L = 8. These results show that the all- 
out spin structure with uniform scalar chirality is indeed stable 
for J/t < 4, while the ferromagnetic phase is stabilized for 
J/t > 5. Our MC simulation does not converge at these very 
low temperatures in the interval 4 < J/t < 6. We traced the 
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FIG. 2: (color online) (a) Energy per site as a function of J/t for 
the all-out, 2in-2out, and ferromagnetic spin structures at p — 1/4 
for L = 512 (solid line) and L = 8 (dashed line). Insets show an 
enlarged view for small J/t. (b) Coupling J/t dependence of Xo and 
5(0) at quarter filling by Monte-Carlo method with L = 8 and very 
low temperature T/t = 1CP 4 . Four mean values of Xo and 5(0) are 
shown for each value of J/t to indicate the lack of convergence the 
MC simulation in the interval the phase transition, 4 < J/t < 6. 
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FIG. 3: (color online) Results of Monte-Carlo simulations at quarter 
filling and J/t = 2. (a) Temperature dependence of specific heat c 
and Xo ■ (b) Snapshot of a spin configuration. Each color corresponds 
to each of the four different sublattices (L = 12, T/t = 0.005). (c) 
Internal energy distribution near the critical temperature for L = 12. 



lack of convergence to an intervening spiral phase, which in 
the L = 8 systems is slightly lower in energy than the fully po- 
larized and the all-out phases for 4.9 < J/t < 6.3. However, 
for larger systems (L > 64), the spiral state becomes unstable. 
Therefore, we conclude that the non-convergent window be- 
tween the all-out and the fully polarized phases is only a finite 
size effect, and there is a direct transition from the all-out to 
the ferromagnetic phase as a function of J 12311 . We also find 
that S(0) is finite while Xo 1S almost zero below J/t ~ 0.5. 
This is a clear consequence of the finite size effect that was 
already discussed in our variational calculation [see Fig|2a)]. 
The intermediate coupling range, 1 < J/t < 4, appears to be 
stable against size effects. 

We now turn to the question of the finite temperature stabil- 
ity of the chiral magnetic phase. Our main results are pre- 
sented in Fig [3] To control finite size effects, we studied 
several system sizes L at fixed coupling strength J/t = 2. 
This value of J/t gives a robust all-out chiral ordering at 
T = 0. As expected from the simple Ising argument in 2D, 
this chiral ordering should persist for a finite range of tem- 
peratures. Starting from the maximum possible value, the 
scalar chirality decreases as a function of temperature, van- 
ishing at T c ~ 0.026t. The snapshot of a spin configuration at 
T = 0.005t in Fig. [3jb) illustrates how the finite-temperature 
fluctuations destabilize the chiral ordering. The calculated 
specific heat curve exhibits a very sharp and symmetric peak 
at T c , which is an indication of a first order phase transition. 
To test this possibility, we analyzed the temperature depen- 
dence of the internal energy probability distribution function, 
Fig. Oc). We found that it has a bimodal distribution for T 
near T c and L = 12. In combination with the specific heat 
behavior, this result indicates that the corresponding thermo- 
dynamic phase transition between the paramagnetic and the 
chiral all-out state is of the first order. 



Magnetic field, H, is an important parameter that can be 
used to control the chiral magnetic states. It couples to both 
local and itinerant electron spins, as well as to the orbital mo- 
tion of electrons, and hence can lead to rich variety of phases. 
For simplicity, here we consider the minimal - Zeeman - 
coupling to the local moments 



Eh 



(4) 



with H = Hz, and neglect the Zeeman and orbital coupling 
to the conduction electrons. This approximation is rigorously 
justified when the Hund's coupling is much smaller than the 
itinerant electron bandwidth. In this case, the electron Pauli 
susceptibility can be neglected relative to the local moment 
susceptibility, while the orbital effect of Berry curvature in- 
duced by the magnetic texture dominates the orbital effect of 
the applied field. 

The T = variational approach for J/t = 2 reveals that 
the all-out structure is distorted by the field in the way that 
interpolates between the all-out and three-in one-out struc- 
ture [Fig. Hlb)]; until it saturates at H sat ~ 0.5t as shown 
in Fig. |Sb). The spins of one of the four sublattices are al- 
ways aligned with the field direction. The gap in the electron 
spectrum closes at the critical field H c /t ~ 0.15 leading to an 
insulator-to-metal transition. Correspondingly, a xy , starts to 
decrease continuously at H = H c , from its quantized value 
for H < H c to zero for H > H sat . In this way, the appli- 
cation of an external field suppresses SQHE in the insulating 
state by inducing a metallic phase with anomalous Hall effect. 

We performed MC simulations in the low field region, H = 
0.025t and H = 0.05i, for L = 8 and L = 12 (for larger 
values of H the variational approach revealed spurious finite 
size effects in systems with L < 64). The results are shown in 
Fig. [5] The position of the specific heat peak (Fig. EJa)) and 
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FIG. 4: (color online) Variational calculation of the ground state of 
% + H z by minimizing the energy over all possible four-sublattice 
spin configurations (L — 2048): a) Hall conductance, a xy , and uni- 
form magnetization along the field direction, m z , as a function of 
the applied magnetic field H/t. (b) Field dependence of the four- 
sublattice spin structure of the variational ground state. 



the onset of the chiral order parameter (Fig. |5ja)) indicate that 
the transition temperature to the chiral phase is suppressed by 
the presence of the magnetic field. The snapshots of the spin 
configurations in each sublattice (identified with color) show 
a result that is consistent with the variational approach: the 
spins are canted towards the three-in one-out configuration, 
with one of the spin sublattices tending to be aligned with the 
applied field (see Figs. |5jb) and (c)). The fluctuations around 
this preferred orientation decrease with increasing H. 
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FIG. 5: (color online) Results of Monte-Carlo simulation at quarter 
filling and J/t — 2, for two values of applied magnetic field, H/t = 
0.025 and 0.05, (a) temperature dependence of specific heat c and 
Xo , (b,c) configurations of all spins classified according to sublattice 
(L = 12, T/t = 0.003), with (b) H/t = 0.025, (c) H/t = 0.05. 

In summary, we have shown that a chiral spin liquid is 
present in the quarter-filled KLM with classical local mo- 
ments on a triangular lattice. This liquid exhibits spontaneous 
quantum Hall effect, which can be tuned by an external mag- 
netic field. The phase is stable with respect to thermal fluc- 



tuations and exists in a wide range of the Hund's coupling 
strengths. It can therefore be realized, e.g., in manganese- 
based materials [24J. It may also be relevant to Na^CoC^ 
|lil|25|] in the limit of quantum local moment, S = 1/2. 
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